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Cretaceous represents one of the hottest greenhouse periods in the Earth's history, but some recent studies 
suggest that small ice caps might be present in non-polar regions during certain periods in the Early 
Cretaceous. Here we report extremely negative 5 ls O values of — 18.12%o to — 13.19%o for early Aptian 
hydrothermal zircon from an A-type granite at Baerzhe in northeastern China. Given that A-type granite is 
anhydrous and that magmatic zircon of the Baerzhe granite has 5 18 0 value close to mantle values, the 
extremely negative 5 18 0 values for hydrothermal zircon are attributed to addition of meteoric water with 
extremely low 5 ls O, mostly likely transported by glaciers. Considering the paleoaltitude of the region, 
continental glaciation is suggested to occur in the early Aptian, indicating much larger temperature 
fluctuations than previously thought during the supergreenhouse Cretaceous. This may have impact on the 
evolution of major organism in the Jehol Group during this period. 



Although the Cretaceous is reputed to have had a hot, greenhouse climate 1 " 4 , its entire history is not well 
documented. Most published isotopic data, paleobiogeography of terrestrial and marine organisms, leaf 
physiognomy, and the distribution of climatic-sensitive sediments have been interpreted to indicate that 
temperatures in the Early Cretaceous were much warmer than present, especially at high latitudes 2,3 . However, 
there are some paleontological records implying a cool climate at that time 5 " 7 . This interpretation is also sup- 
ported by other indications of cool intervals from the late Barremian to the early Albian of the Early Cretaceous 6,7 , 
with an estimated mean air temperatures comparable to the present day "icehouse" — about 10 ± 4°C at mid- 
latitudes (~42°N) in Asia 8 , implying small polar ice caps. Concrete evidence for these ice caps is sparse, however, 
because the subsequent oscillation between greenhouse and icehouse conditions may efficiently eliminated the 
surface records of paleo-glaciation 5 . Here we provide new 5 ls O data of zircon from the Baerzhe A-type granite in 
northeastern China at the paleolatitude of ~ 45 °N, which indicate a contribution from continental glacial 
meltwater (Fig. 1). 

The Baerzhe A-type granite is located in southeastern Great Khingan Mountains, northeastern China (Fig. 1), 
which intruded into the Late Jurassic and Early Cretaceous igneous rocks, in the eastern part of the Central Asian 
Orogenic Belt 9 . The Late Jurassic volcanic rocks consist of andesite, andesitic tuff, and intermediate- acid volcanic 
clastic rock 10 . The Early Cretaceous igneous rocks include rhyolites and andesitic tuffs that unconformably 
overlying the Late Jurassic volcanic rocks. The giant Baerzhe Zr-REE-Nb deposit is hosted in the albite-rich 
phase, which occurs in the upper section of the Baerzhe granite, while the lower section is barren (or weakly 
mineralized). Pegmatites were discovered on the top of the mineralized granite 11,12 . The granite consists mainly of 
microcline, quartz, arfVedsonite and albite, whereas the proportions, grain sizes and crystal forms of rock- forming 
minerals differ in the mineralized and barren counterparts. The mineralized granite contains much more aegirine 
and albite and less arfVedsonite, moreover, the crystals are larger and slightly more euhedral, which are in striking 
contrast to the barren granite. The occurrence of abundant miarolitic cavities suggests a shallow level emplace- 
ment (< 2 km) 13 and highly evolved 9,14 . Zircon is an important ore mineral in the Zr-REE-Nb deposit, with both 
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Figure 1 | Geographical positions of the Baerzhe intrusion and the deposits of Jehol Biota. Inset is the distribution of the Jehol ecosystem 
(radial circle) in the paleogeographic map of eastern Asia in the Early Cretaceous, modified from Zhou et al 57 . The Jehol Group is comprised of the Yixian 
and Jiufotang Formations, which crop out in western Liaoning, northern Hebei and southeastern Inner Mongolia 57 . The paleolatitude of the Baerzhe 
A-type granite is ca. 45°N. Abbreviations refer to major tectonic divisions: EUR, Europe; INC, Indo-China; IND, India; JUN, Junggar; KAZ, Kazakhstan; 
LH, Lhasa; MON, Mongolian; NCB, north China; QA, Qiadam; QI, Qiangtang; SCB, south China; SH, Shan Thai; SIB, Siberian; and TAR, Tarim. 



magmatic and hydro thermal origins. The hydrothermal zircon we 
identified in Baerzhe is characterized by features originated from the 
magmatic-hydrothermal transition systems of the highly evolved 
granitic pluton 9 , which include: (1) tetragonal dipyramidal in mor- 
phology that is a relatively low- temperature zircon type from evolved 
alkaline granite 1516 ; (2) murky and featureless in textures and LREE- 
rich and high common lead in compositions are common character- 
istics of zircon grains that directly crystallized from Zr- saturated 
hydrothermal fluids 17,18 . 

Results 

Magmatic zircon grains from the A-type granite occur as crystals of 
prismatic morphology with oscillatory zonation in cathodolumines- 
cence (CL) images and are depleted in LREE (Supplementary Table 
SI). They have concordant 206 Pb/ 238 U ages of 123.9 ± 1.2 Ma (2a, 
MSWD = 0.65, n = 17) (Fig. 2; Supplementary Table S2). They have 
positive 5 ls O values of 2.79 to 5.10%o (Fig. 3; Supplementary Table 
S3), slightly lower than values of 5.3 ± 0.3%o for normal mantle 
zircon 5 ls O values 19 . They have e H f{t) values ranging from 1.51 to 
12.6 (Supplementary Table S4). Quartz phenocrysts from the granite 
exhibit an average 5 ls O value of 5.37 ± 0.54%o (la) (Supplemen- 
tary Table S5). When paired with a weighted average 5 ls O value of 
3.47 ± 0.1 l%o (la) for the magmatic zircon, this yields an O isotope 
temperature of 1030°C (Supplementary Fig. S8) 20 . This temperature 
corresponds to the highest value for homogenization temperatures of 
primary melt inclusions (750 to 1030°C) in quartz phenocrysts 12 , 
suggesting O isotope equilibrium between the magmatic quartz 
and zircon in the granite. 



Hydrothermal zircon grains from the Baerzhe ore deposit occur as 
crystals of tetragonal dipyramids, which are murky in CL images and 
enriched in LREE (Supplementary Table SI). They have concordant 
206 pb/ 238Tj age 0 f 12 3.5 ± 3.2 Ma (2a, MSWD = 0.37, n = 7) (Fig. 2; 
Supplementary Table S2), and have s H f(t) values ranging from 2.54 to 
7.62 (Supplementary Table S4). Most analyses for the hydrothermal 
zircons yield extremely negative 5 ls O values of — 18.12%o to 
— 13.19%o (Fig. 3; Supplementary Table S3), which are significantly 
lower than magmatic zircon from other Cretaceous granite in the 
adjacent area 21,22 . To our knowledge, — 18.12%o is the lowest 5 ls O 
value of zircon from granite so far reported 23 " 28 . It is significantly 
lower than negative 5 ls O values of — 14.4 to — 10.0%o for magmatic 
garnet from granite 27 and —10.9 to — 7.8%o for metamorphic and 
magmatic zircons from metabasite and metagranite 25 ' 28 in the Dabie- 
Sulu orogenic belt. Such negative 5 ls O values from the Dabie-Sulu 
rocks have been attributed to high-temperature hydrothermal altera- 
tion in a rift tectonic zone by continental glacial meltwater prior to 
the Snowball Earth event 26 ' 27 . Since O diffuses very slowly in zir- 
con 29 ' 30 , crystalline zircon is capable of preserving the original 5 ls O 
signature of its source 23 ' 25 . The extremely negative ls O values for the 
hydrothermal zircon from the highly evolved Baerzhe pluton indi- 
cate that unusual ls O- depleted water was incorporated into the Zr- 
saturated hydrothermal fluid in which the hydrothermal zircon was 
crystallized. 

Discussion 

Fluid 5 ls O values responsible for the negative 5 ls O hydrothermal 
zircon grains are estimated to be lower than — 15.23%o at 600°C, 
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Figure 2 | Zircon U-Pb ages for the Baerzhe A-type granite and its coincidence with the Jehol Biota deposits, drop of global sea-level and negative shift 
of carbon isotopes for Pacific sediments. The ages range from 124.1 ± 0.3 Ma to 122.9 ± 0.7 Ma of a tuff in the main fossil bed of Jehol Biota in the Yixian 
Formation is highlighted with gray 58 . The low value on the eustatic curve in the early Aptian is attributed to glacioeustasy 37 . The carbon isotopes curves 
also show a significant negative excursion at ca. 124 Ma 38 . 



according to the O isotope fractionation factor between zircon and 
water 20 . If hydrothermal zircon is crystallized in hydrothermal fluid 
consisting of a mixture between a magmatic exsolved fluid and a 
surface water, 5 ls O value for the participated surface water should 
be lower than the lowest measured 5 ls O value of zircon that we 
report. We estimate that the water that caused the extremely ls O 
depletion of hydrothermal zircon could have been as low as — 17%o 
to — 34%o, assuming 10% to 50% of the fluid was magmatic exsolved 
fluid. The only significant terrestrial reservoir with such negative 
5 ls O values is the highly fractionated meteoric water and its frozen 
product (continental ice sheet) 26,31 . 

The most plausible source of such negative 5 ls O values can be 
expected from continental glacier melt water. The Baerzhe A-type 
granite is located in NE Asia at mid-latitude and not far away from 
the Pacific Ocean (< 700 km) since the Jurassic (Fig. 1) based on the 
paleomagnetic results 32 . The current altitude of the Baerzhe region is 
850 m. This region did not experience major denudation as indicated 
by Late Jurassic calderas remnants 33 , such that the paleo- altitude of 
the Baerzhe region is inferred to be ~ 2 km or less 34 , which was 
exclude of intra-continental high mountains. 5 ls O values for pre- 
sent-day precipitation at the Baerzhe region is estimated to be about 
— 10%o 35 , and it is likely to be more or less the same in the Early 
Cretaceous. Such low 5 ls O meteoric water (< — 17%o) is seen in high- 
latitude regions, i.e., modern meteoric precipitation of the Greenland 
ranges from — 15%o to — 35%o 35 . These kinds of regions usually have 
very low precipitation. Given that Baerzhe is located several thou- 
sand kilometers away from such low 5 ls O region, with higher pre- 
cipitation and much higher 5 ls O, it is not likely that water from high 
latitudes can be transported to Baerzhe through rivers and keeps the 
very low 5 ls O values. Glacier is the only plausible way that transports 
water for thousands of kilometers and keep the extremely low 5 ls O 
values. Therefore, continental glacial meltwater is the most likely 
source for the negative 5 ls O water ( — 17%o to — 34%o), responsible 



for the ls O-depleted Baerzhe hydrothermal zircon at mid-latitudes 27 . 
The occurrence of glacier meltwater incorporated in the hydro- 
thermal zircon can be interpreted to imply mean annual tempera- 
tures of < — 7.5°C in the target region (Fig. 4). 
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Figure 3 | Plot of £nf(t) versus 5 18 0 values for magmatic and 
hydrothermal zircon in this study, showing remarkable shifts of oxygen 
isotopes corresponding to incorporation of continental glacial meltwater 
into the alkaline magma. VSMOW denotes the Vienna Standard Mean 
Ocean Water. The normal mantle zircon 5 ls O value is 5.3 ± 0.3%o 19 . The 
consistent s H f(t) values and great shift of 5 ls O values from the magmatic 
zircon to the hydrothermal zircon were caused by the incorporation of 
continental glacial water. Note that the three middle spots of magmatic 
zircon show a minor hydrothermal influence. 
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Figure 4 | The mean annual temperature (MAT, °C) can be estimated 
from the 5 ls O values (%o, VSMOW) of hydrothermal fluid that yielded 
the hydrothermal zircon. If the hydrothermal fluid was composed of 
different proportions of surface water (from 90% to 10%) and magmatic 
water, the estimated mean annual temperature of lower than — 7.5 °C, 
using the following equation: 5 ls O = 0.64 X T (°C) - 12.8 56 . 

In principle, glacial meltwater may be incorporated into hydro - 
thermal zircon either by direct mixing with magmatic fluid, or by 
assimilation of supracrustal rocks that have previously been altered 
by the glacial meltwater at high temperatures. The latter requires 
interaction between glacial meltwater and granitic magma during 
the emplacement, imparting the extremely negative 5 ls O values to 
cooling magmas to generate the similar 5 ls O values of hydrothermal 
fluid. Although the O isotopes of hydrothermal zircon are controlled 
by both magmatic water and glacial meltwater, the budget for water- 
insoluble trace elements in zircon is dominated by the granitic melts. 
In contrast to the dramatic difference in O isotopes between the 
hydrothermal and magmatic zircon, Hf isotopes of hydrothermal 
zircon fall within the range of those of magmatic zircon (Fig. 3). 
This suggests that the hydrothermal zircon gained extremely nega- 
tive 8 18 0 without significant changes in Hf isotope composition, 
which does not favor the assimilation of supracrustal rocks during 
magma emplacement. The whole rock Hf-Nd isotopes analyses also 
suggest that the Baerzhe granitic magma was derived from partial 
melting of juvenile crust, without significantly supracrustal assimila- 
tion (Supplementary Fig. S5). Therefore, extremely low 5 ls O water 
participated in shallow intrusions 24 ' 36 is the favorable explanation for 
Baerzhe hydrothermal zircon. 

The occurrence of continental glacial meltwater indicates the 
presence of an exceptionally cold paleoclimate locally in north- 
eastern Asia in the early Aptian. This is supported by global sea- 
level changes, carbon isotopes in oceanic sediments and some 
paleontological data from Jehol Biota. Sea-level records from 100 
to 155 Ma suggest a warm global climate during the Earliest 
Cretaceous, followed by dramatic sea-level falls, which was attrib- 
uted to glacioeustasy from the late Barremian (ca. 128.3 Ma) to the 
early Aptian (ca. 123.3 Ma) 6 ' 37 . In addition, there is a major nega- 
tive shift of carbon isotopes for the Pacific carbonate section at 
~ 124 Ma 38 , indicating abrupt methane emission that may have 
been triggered by oceanic deglaciation 39 . All these suggest that the 
continental ice sheet may have developed locally in northeastern 
Asia, and supplied deglacial meltwater to mid-paleolatitude regions 
(~45°N) in the early Aptian. 



Consistently, cool paleoclimate is supported by paleontological evid- 
ence, including plant fossils such as fossil wood genus Xenoxylon 40 , 
insect fossils such as stoneflies and raphidiopterans 41 , and the absence 
of thermophilic reptiles such as crocodilians 8 . Most significantly, the 
recent discovery of a gigantic feathered theropod dinosaur with long 
feathers in the Early Cretaceous from Jehol Biota was proposed to be 
analogous to some large mammal taxa, e.g. mammuthus primigenius, 
in ice age, in terms of developing long integumentary coverings as an 
adaptation to an unusually cold environment 42 . 

Previous studies demonstrate the presence of cold intervals during 
the greenhouse Mesozoic 5,6,37 , but the documentation of continental 
glaciation in mid-latitude regions is somewhat unexpected. The Jehol 
Biota lasted for about 10 million years, whereas it is likely that this 
early Aptian ice age represents only a short-lived cold interval during 
existence of the Jehol Biota. Our study thus encourages more high- 
resolution geochemical and other studies which may reveal more 
cold intervals during the Mesozoic greenhouse era. The exceptionally 
cold climate in the Early Cretaceous in northeastern Asia is also 
inconsistent with abundant evidence supporting a much warmer 
climate in low-latitude regions 1 ' 4 ' 5,43 . This can be explained by the 
hypothesis that the climate regime of the Earth during the 
Mesozoic is different from the present, i.e., a much steeper pole-to- 
equator temperature gradient in the Mesozoic than today 5 . 

The increasing evidence supporting large temperature fluctua- 
tions during the greenhouse Cretaceous has implications for our 
understanding of the evolution of the Mesozoic ecosystem. There 
is some evidence supporting the wide distribution of feathery or fury 
coverings in dinosaurs, pterosaurs, and mammals, the three dom- 
inant Mesozoic groups 44 ' 45 . The relative development of these integu- 
mentary coverings is possible to be linked with the tremendous 
temperature fluctuations, though the available data is far from 
enough yet to make such an inference. 

In conclusion, the extremely low 5 ls O values for the hydrothermal 
zircon presented in this study favor an explanation by the incorp- 
oration of continental glacial meltwater rather than intracontinental 
water or the assimilation of previously altered rocks. The required 
glacial meltwater (— 17%o to — 34%o) was corresponding to mean 
annual temperature of lower than — 7.5°C, suggesting that contin- 
ental ice sheet would had developed to the mid-latitude and low- 
altitude regions in NE Asia in the Early Cretaceous. This is a surpris- 
ing result, because the Cretaceous is known to be one of the hottest 
periods in the Earth's history. More broadly, such extremely cold 
climates in the supergreenhouse Early Cretaceous might have 
brought significantly forces on the evolution of the Mesozoic ecosys- 
tem, including the Jehol Biota. 

Methods 

Zircon grains were separated using the standard density and magnetic separation 
techniques. Both hydrothermal and magmatic zircon grains were selected from 
mineralized and barren Baerzhe alkaline granites. These grains were handpicked 
under a binocular microscope and mounted in an epoxy resin disc, then polished to 
near half section to expose internal structures. Transmitted and reflected-light 
microscopy, CL and scanning electronic microscope (SEM) investigations were 
carried out before in-situ U-Pb dating and O-Hf isotopic analyses on inclusion-free 
domains. 

SIMS U-Pb dating. U-Pb dating was conducted using a Cameca secondary ion mass 
spectrometer (SIMS) 1280 at the Institute of Geology and Geophysics, Chinese 
Academy of Science (IGGCAS) (Supplementary Table S2). Analytical procedures are 
the same as those described by Li et al. 46 . The 0 2 ~ primary ion beam was accelerated at 
13 kV, with an intensity of ca. 8 nA. The ellipsoidal spot is about 20 X 30 um in size. 
Positive secondary ions were extracted with a 10 kV potential. The Pb/U calibration 
was performed using TEM zircon standard with an age of 417 Ma 47 . Data processing 
was carried out using the Isoplot 4.11 programs of Ludwig 48 , and the 204 Pb-based 
method of common Pb correction was applied. The uncertainties in single analyses 
are cited as la, and the weighted mean ages are quoted at the 95% confidence level 
(2a). 

SIMS O isotopes analyses. The oxygen isotopic compositions of zircon were 
measured using the same Cameca IMS- 1280 SIMS at the IGGCAS (Supplementary 
Table S3), with analytical procedures similar to those reported by Li et al. 49 . The Cs + 
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primary ion beam was accelerated at 10 kV, with an intensity of ca. 2 nA (Gaussian 
mode with a primary beam aperture of 200 um to reduce aberrations) and rastered 
over a 10 urn area. The spot is about 20 um in diameter (10 um beam diameter + 
10 um raster). The normal incidence electron flood gun was used to compensate for 
sample charging. Negative secondary ions were extracted with a — 10 kV potential. 
Oxygen isotopes were measured in multi-collector mode using two off-axis Faraday 
cups. The mass resolution used to measure oxygen isotopes was ca. 2500. The nuclear 
magnetic resonance probe was used for magnetic field control with stability better 
than 2.5 ppm over 16 hours on mass 17. The instrumental mass fractionation factor 
(IMF) is corrected using zircon 91500 standard with a 5 ls O value of 9.86 ± 0.1 l%o 50 . 
Penglai zircon megacrysts with a reference 5 ls O value of 5.31 ± 0.10%o were used for 
further external adjustments 51 (Supplementary Fig. S3). Measured 18 0/ 16 0 ratios were 
normalized by using Vienna Standard Mean Ocean Water compositions (VSMOW, 
18Q/16Q _ 0.0020052), and then corrected for the instrumental mass fractionation 
factor (IMF) as follows: (§ 18 0) M = [( 18 O/ 16 O) M /0.0020052-l] X 1000 (%„), IMF = 

(0 18 0) M (standard) " (o 18 0) VS MOW> 0 18 0 Samp l e = (§ 18 0) M + IMF. 

LA-MC-ICPMS Hf isotopes analyses. Hf-isotope analyses were carried out in-situ 
using a Nu Plasma HR MC-ICP-MS with a GeoLas 2005 excimer ArF laser- ablation 
system, at the State Key Laboratory of Continental Dynamics, Northwest University, 
Xi'an, China. The Nu Plasma HR MC-ICP-MS is a second- generation double- 
focusing MC-ICP-MS with three ion counters and twelve faraday cups. The GeoLas 
2005 laser-ablation system consists of a COMPexPro 102 ArF excimer laser 
(wavelength of 193 nm, maximum energy of 200 mj, and maximum pulse rate of 
20 Hz) and a GeoLas 2005 PLUS package. Lu-Hf isotopic analyses were obtained on 
the same zircon grains that were previously analyzed for O isotopes, with ablation pits 
of 44 um in diameter, ablation time of 26 seconds, repetition rate of 7 Hz, and laser 
beam energy density of 15 J/cm 2 . The analytical procedures were similar to those 
described by Yuan et al. 52 . Zircon 91500, Monastery and GJ-1 were used as external 
standards in this analysis, and our 176 Hf/ 177 Hf ratios for the three zircon standards 
(Supplementary Fig. S4) are in good agreement with reported values 52 . Hf isotopic 
data for the hydrothermal and magmatic zircon are listed in Supplementary Table S4. 

Estimation of mean annual temperatures. Oxygen isotope ratios provide a powerful 
tool for understanding the thermal, magmatic and fluid history 53 ' 54 . Yuan and Zhang 55 
reported 5 ls O values for whole rock and quartz, K-feldspar, zircon and albite from the 
Baerzhe granite (Table S3, Fig. S6). The highly negative 5 ls O values for whole-rock, K- 
feldspar and zircon suggest O isotopic disequilibrium between minerals in the granite. 
The huge A 18 OQ uartz _ K .f eldspar values show that meteoric water has interacted with the 
granite (Supplementary Fig. S7). 

The weighted average value of the 34 5 ls O analyses for the magmatic zircon 
(excluding three mixed spots) is 3.47%o ± 0.49%o (la), which is significantly lower 
than the 5 ls O range of 5.3%o ± 0.3%o for zircon in high-temperature equilibrium with 
the normal mantle 19 . The statistical mean value of the 10 quartz 5 ls O analyses is 5.37 
± 0.54% 0 (la). According to the theoretical calculations of O isotope fractionation 20 , 
equilibrium O isotope fractionations between quartz and magmatic zircon are 
1.90 ± 0.73%o (la) at magmatic temperatures. From the following fractionation 
equation 20 : 

10 3 lna quartz . zircon = 0.72 x 10 6 /T 2 +4.26 x 10 3 /T- 1.79 

At equilibrium 10Tna quartz „ zircon » A 18 O quartz . zircon » 5 18 O quartz - § 18 O zircon , we can 
estimate an equilibrium temperature (T) of ca. 1030°C (Supplementary Fig. S8). This 
temperature is only slightly higher than the homogenization temperature of 750- 
1030°C from primary melt inclusions in quartz phenocrysts 12 . 

The minimum 5 ls O value of 31 hydrothermal zircon grains is — 18.12 ± 1.5%o, 
which, from the theoretical calculations of zircon-water fractionation 20 , indicates 
hydrothermal fluid with a maximum 5 ls O value of — 15.23%o. The equilibrium 
temperature was assumed at 600°C, which consistent with the homogenization 
temperatures of 475 to 650°C from melt-fluid inclusions in the pegmatite of the 
Baerzhe deposit 12 . 

10 3 lna zircon _ wate r = 3.76x 10 6 /T 2 -9.03 x 10 3 /T + 2.52 

At equilibrium 10 3 lna zircon . water ~ A 18 O zircon . water ~ 5 18 O zircon - 5 18 O water . If T = 
600°C, 10 3 lna zircon _ water = -2.89%o. Using these values, we obtain the 5 ls O value of 
— 15.23%o for the hydrothermal fluid. Whereas, the hydrothermal fluid that yielded 
the hydrothermal zircons should be composed of surface water (SW) and magmatic 
water (MW). The A-type granite normally is anhydrous, which means that only a 
limited amount of magmatic fluid was exsolved from the Baerzhe magma. Assuming 
that the proportion of incorporated surface water (w\) is < 90%, and the primary 
magmatic water is (w 2 ) is > 10%, the § ls O value of the hydrothermal fluid is calcu- 
lated as: 

(5 18 O hf = wi x (5 18 O sw + w 2 x 3 l8 0 M w 

where, w 1 + w 2 = 1. 

The relationship between 5 ls O value of continental surface water and annual 
temperature (T, °C) 56 is as follows: 

(5 18 O sw = 0.64 xT- 12.8 

Furthermore, the magmatic fluid exsolved from the Baerzhe magma should be in O 
isotopic equilibrium with the magmatic zircon. From the theoretical calculations of 



zircon-water 20 , we obtain a § ls O M w value of 6.36 ± 0.52% for the magmatic water at 
600°C. 

10 3 lna zircon _ water = 3.76 x 10 6 /T 2 -9.03 x 10 3 /T + 2.52 

~<5 18 Omagmatic z i rcon — ^OjMW 

= -2.89 

We can now estimate the mean annual temperature using the 5 ls O value of the 
hydrothermal fluid, as show in Fig. 4. 

It is possible that the incorporation of surface water at the magmatic-hydrothermal 
transition (ca. 600° C), varied in different parts of the pluton. Assuming the lowest 
measured 5 ls O value of — 18.12 ± 1.5%o in the hydrothermal zircon is composed of 
90% surface water in the hydrothermal fluid, the corresponding mean annual tem- 
perature is ca. — 7.5°C. The range of 5 ls O values of the hydrothermal zircon (—18.12 
± 0.15%o to — 13.19 ± 0.14%o) could be due to variably amounts of surface water 
(continental glacial meltwater). 
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